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Utilizing the Langevin equation for the linear u-model we 
investigate the effect of friction and white noise on the evolu- 
tion and stability of collective pionic fields in energetic heavy 
^ — „ion collisions. We find that the smaller the volume, the more 
QV^table transverse (pionic) fluctuations become on a homogeneous 
(3\lisoriented chiral field background (the average transverse mass 
r- H(m,() increases). On the other hand the variance of increases 
Ljeyen more, so for a system thermalized in an initial volume of 10 
Q£m 3 about 96% of the individual trajectories enter into unsta- 
■^^ble regions (to? < 0) for a while during a rapid one dimensional 
expansion (to = 1/m/c). In contrast the ensemble averaged so- 
^ "lution in this case remains stable. This result supports the idea 
of looking for disoriented chiral condensate (DCC) formation in 
^individual events. 
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fN) Ultrarelativistic heavy ion collisions offer the possibility 
■^"t<j> study hadronic matter at high initial temperatures and 
•Oenergy densities. Because of their small masses pions con- 
£^tribute most dominantly to the large multiplicity of pro- 
educed particles. It was speculated that some of the observed 
^^pions are produced by the coherent decay of a semiclassi- 
pion field The idea of so called 'disoriented chi- 

Qfj fl? condensates' (DCC) was then made widely known due 
QJto Bjorken, Kowalski and Taylor ||. They speculated that 
^^events can occur in which the classical pion field is oriented 
^.aiong a single direction in isospin space. 

k>( | Rajagopal and Wilczek Q have proposed that in a rapid 
5_jquench through the second order chiral QCD phase transi- 
tion such disordered chiral configurations may emerge. Un- 
stable, exponentially increasing long wavelength pion fields 
can develop during the 'roll down' period of the order param- 
eter according to the pure classical equation of motion j4||| . 
The spontaneous growth and subsequent decay of these con- 
figurations would give rise to large collective fluctuations in 
the number of produced neutral pions compared to charged 
pions, and thus could provide a mechanism explaining a fam- 
ily of peculiar cosmic ray events, the Centauros |6) . A deeper 
reason for these strong fluctuations lies in the fact that all 
pions are assumed to sit in the same momentum state and 
the overall wavefunction can carry no isospin [Q . 

However, the proposed quench scenario assumes a highly 
non-equilibrium initial configuration for the phase transition 
to occur, namely that the potential governing the evolution 
of the order parameter and the long wavelength modes im- 
mediately turns to the classical one governing the vacuum 
structure at zero temperature. This represents a very dras- 
tic assumption as the soft and classical modes completely 
decouple in an ad hoc manner from the residual thermal 



fluctuations at or slightly below the critical temperature. It 
is not likely to happen in an ultrarelativistic heavy ion col- 
lision. 

An alternative scenario was suggested by Gavin and 
Muller H who proposed to use instead of the bare potential 
the effective one-loop potential including the thermal fluc- 
tuations on the mean-field level. If the system cools rapidly 
enough (due to longitudinal (D=l) or radial (D=3) expan- 
sion) and, similar to the quench scenario, if the order pa- 
rameter is small enough at the onset of the evolution below 
the critical temperature, DCC might emerge as well. The 
latter assumption, however, has also been criticized: if the 
order parameter stays in thermal contact with the fluctua- 
tions giving rise to the effective 1-loop potential, then one 
also has to allow for the thermal fluctuations in the initial 
conditions p|J[o[|. Quenched initial conditions seem to be 
statistically unlikely. 

In this letter we address the question of the 'likcliness' 
of an instability leading potentially to a DCC event in the 
light of different initial conditions as raised by Blaizot and 
Krzywicki Jjj|| . 

In a simplified model we study the influence of the ther- 
mal modes ('fluctuations') on the off-equilibrium evolution 
of the order parameters. These chiral fields are treated semi- 
classically within the chiral 0(4) cr-model. The long wave- 
length modes, especially the zero modes used as order pa- 
rameters, represent an open system which constantly inter- 
acts with the thermal fluctuations (the 'reservoir'). On the 
one-loop level the interaction with the hard modes will gen- 



erate an effective mass term \\T 2 [p-|l0|, leading to the 



effective Hartree-Fock potential. It was recently shown in 
detail that in the (/) 4 -theory hard modes can be integrated 
out on the two-loop level leading to dissipation and noise in 
the quasi-classical limit for the propagation of the long wave- 
length fields The resulting equations of motion are of 
Langevin-type. In the weak coupling limit, the friction coef- 
ficient r\ is directly related to the on-shell plasmon damping 
rate, r\ = 27 P ;. The noise term shows up as an imaginary 
part in the effective action and is related to the friction via 
the fluctuation-dissipation theorem. The interplay of noise 
and dissipation guarantees that the soft modes eventually 
become thermally populated on the average. 

We propose to study the following Langevin equations of 
motion for the order parameters $ a = y J d 3 x (j) a (x,t) in a 
volume V 



1 



(1) 



with $ a = (cr, 7T , 7r 2 , 7r 3 ) being the chiral meson fields and 



i2 = a(*8+5>? 



(2) 



Here r is the proper time of the expanding system and the 
'dot' denotes the derivative with respect to r. The above 
equations assume a D-dimensional scaling expansion [|s| — p~of| . 

Before presenting our results some comments are in or- 
der: We use the standard parameters = 93 MeV for the 
pion decay constant, to w = 140 MeV for the pion mass, and 
A = 20 for the coupling constant. With this choice we are 
obviously in the strong coupling regime, so our conclusions 
drawn from the investigation of the weak coupling regime 
]l2|] are rather a motivation for our present usage of simple 
effective terms for noise and dissipation. 

We treat the dissipation term as Markovian, which as- 
sumes a clear separation among the timescales of the hard 
and soft modes. In the semiclassical Markovian approxima- 
tion the noise is effectively white and at two-loop level it is 
Gaussian lOl, 
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The phenomenon of long wavelength DCC amplification 
occurs in periods when the transverse mass squared m 2 be- 
comes negative. This happens during a quench scenario or 
rapid cooling. For the sake of simplicity we treat 77 as a 
constant, obtained using (Q) and m/T « 1 throughout the 
evolution. At T = T c = ^2/ 2 - 2m%/\ = 123 MeV the 
friction 77 = 2^( v i = 2.2 (fm/c) -1 is rather strong. Therefore 
we also investigate scenarios with 1/2 and 1/4 of this value. 

Aside from a theoretical justification one can regard the 
Langevin equation as a practical tool to study the effect of 
thermalization on a subsystem, to sample a large set of pos- 
sible trajectories in the evolution, and to address also the 
question of all thermodynamically possible initial configu- 
rations in a systematic manner. Applying this we are able 
to study the up to now unknown influence of thermal fluc- 
tuations on the growth of disoriented chiral domains. The 
friction coefficient 77 helps to prolong the evolution of the 
order parameter, as it slows down the rolling down of the 
effective potential. On the other hand the noise term leads 
to a subsequent thermalization making DCC formation less 
likely, but allowing also for fluctuations on an event by event 
basis. Thus we need to study how fast the system cools and 
destabilizes and how fast it can thermalize. 

The scaling expansion and cooling of the system is de- 
scribed by the equations 
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where T is the temperature, V the volume and rj the friction 
coefficient. 

Finally we have to specify the friction coefficient 77 for the 
a and pion field. The on-shcll (u> = to) plasmon damping 
rate for standard (/> 4 -theory arising from the 'sunset' diagram 
Jl||.[l2j can be easily taken over to the 0(4)-model, assuming 
that all four masses m for the fluctuations (i.e. hard quanta) 
are equal: 



27 P 
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A 2 — fs P (l 



dtj^j defines the Spence function. 



(4) 



Ad- 



where fs P (x) = - 
mittingly this 'choice' is only a crude estimate as the zero 
modes do not evolve on-shell during the (possibly unstable) 
evolution. Thus the dissipation and noise correlation should 
better be described by nonmarkovian terms including mem- 
ory effects. In addition, the 0(4) transverse and longitudinal 
mass for the fluctuations |§-[Io|, 



= A($ 2 J +5> 2 
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In the framework of the model described so far we in- 
vestigate different evolution scenarios for the order param- 
eters <fr a . For comparison we calculate the pure classical 
(7/ = 0, T = 0) and the one-loop annealing scenario (77 = 0) 
for the k = modes only with quenched initial condition 
($ a = 0, $ a = 0, a = . . . 3). These scenarios contain long 
and highly unstable periods of the evolution during a rapid 
(to = 1/to/c) one-dimesional (D = 1) scaling expansion. 
This can be inspected in the lower part of Fig.l where the 
quantity /it=sign(TO 2 )-\/|TO 2 | is plotted as a function of time. 

In the two-loop motivated Langevin scenario we let the 
system thermalize at temperature T = T c for 10 fm/c from 
the quenched initial condition and then switch on the one- 
dimesional expansion. The outcome of thermalization de- 
pends on the volume occupied by the k — modes - in 
agreement with the equipartition theorem. The middle part 
of Fig.l shows the average evolution of 1000 trajectories, 
each propagated according to the Langevin equation ([!]) , in 
an initial volume of Vb = 1, 10, 100 and 1000 fm 3 . Because 
of the high friction coefficient 77 the equipartition is set after 
the first few fm/c. We find that the smaller the volume, the 
larger to 2 . The volume dependence of the thermalized chi- 
ral order parameter fields can be analyzed using the virial 
theorem 



are not really equal. 



^ 8H . 8H _ 



d<£ a 



(7) 



2 



In the m n = limit analytical results can be obtained as 

3 



= ^ (l ± VI + 16T/AF/2) 
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with / = \{T 2 — T 2 ). The large volume limit of this result 



f + $j + 0(l/V') T<T C 



a=0 



^ + o(i/v a ) 



T = T C (9) 
T >T r 



From this m| = A(£ = (1/^) follows for T < T c 

in accordance with the Goldstone theorem. It is also enlight- 
ening to calculate the energy per chiral degree of freedom 
extra to the mean field energy 
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l + 0(l/V) T>T C 



which reflects the change of the effective number of collective 
chiral degrees of freedom at the phase transition. 
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FIG. 1. ^t(t) for the most unstable event of 1000 with 
V = 10/m 3 , for the average with V = 1, 10, 100 and 1000/m 3 
(from top to bottom), and for the pure and annealing scenarios. 

At the end of the expansion \xt relaxes to the value of 
m n = 140 MeV, as can be seen from the curves in the mid- 
dle part of Fig.l. These ensemble averaged curves do not 



show any significant period of instability. This result is in 
agreement with those of Randrup |l4| ] for a one-dimensional 
expansion from a thermalized initial condition. 

This lack of instability in the averaged evolution does not 
mean, however, that DCC formation cannot be expected in 
heavy ion collisions: Using the Langevin equation we are 
able to explicitely single out particular evolutions which are 
the most unstable. The upper part of Fig.l presents such 
an evolution preheated in a Vq — 10 fm 3 initial volume at 
T = T c . Here, inspite the initial thermalization and ongoing 
noise during the expansion, quite significant unstable peri- 
ods develop. For quantifying the strength of instability we 
define the quantity |3 



G = 



\m t \Q(-m 2 t )dt. 



(11) 



The amplification of small amplitude instabilites with k 
() is then expG. For the particular event shown in Fig.l 
G = 4.7371. The distribution of this quantity P(G), is 
shown in Fig. 2 for an initial volume of Vq = 10 fm 3 , the 
1000 individual trajectories are all preheated to T = T c . 
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FIG. 2. 

V = 10 fm 3 . 



Distribution of the enhancement factor G for 
The most unstable event has G = 4.7371. 



We observe that only in 4.3% of all cases an unstable 
period is missing. The average trajectory shows, however, 
no signal, because the unstable periods occur at different 
times. Therefore the middle part of Fig.l does not allow for 
drawing a conclusion. 

In order to review the tendencies Table 1 summerizes the 
most important properties of the distribution of the amplifi- 
cation factor in the Langevin scenario. The fraction of events 
with no instability at all are written in the second column 
for the various initial volumes. In these cases no DCC signal 
can develop at all. The average growth factor (G) and its 
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maximum within 1000 events, G ma x, 
and fourth columns, respectively. 



are shown in the third 



In the late phase of the expansion the friction (|J) is over- 
estimated by employing m/T = 1, as m/T significantly 
increases due to the cooling and the thermal fluctuations 
decouple. In fact applying a smaller friction the unstable 
oscillating periods appear for a longer time similar to the 
annealing scenario. On the other hand a variation of the 
friction parameter, 77, to its 1/2 or 1/4 value does not change 
significantly the characteristics of the distribution P(G) (cf. 
Table.l). 

In conclusion we have investigated the chiral 0(4) model 
in a Langevin scenario, motivated by the two-loop results of 
the non-equilibrium field theory. In this approach the k = 
order parameter fields are coupled to a thermal bath via 
mass, friction and noise terms. We preheated the system at 
T = T c using its own dynamics and then switched over to a 
one-dimensional scaling expansion. Average and statistical 
properties of individual solutions of the Langevin equations 
were studied with the emphasis on such periods of the time 
evolution when the transverse mass becomes imaginary and 
therefore an exponential growth of unstable fluctations in 
the collective fields can be expected. 

We have found that in different realistic initial volumes 
ranging from 1 to 1000 fm 3 , where the average evolution 
does not show any sensible instability, individual events lead 
to sometimes significant growth of fluctuations. Both the 
strength (G w 2 - 5) and the probability (50 - 90%) of such 
events are remarkably high. The most extreme events are 
quite similar to the predictions of the quenched models. All 
these estimates are optimistic, because the back reaction of 
the soft but / modes will lead to a shortening of the 
unstable periods Jig ], 

Our findings support the idea of looking for DCC forma- 
tion experimentally in individual events. 
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P(G = 0) 


(G) 


Gmax 


1 


% 


3.0983 


7.1163 


10 


4.2 % 


0.8540 


4.7371 


100 


40.3 % 


0.2728 


2.0309 


1000 


40.4 % 


0.1082 


0.6780 


1 


% 


3.0883 


8.5129 


10 


7.0 % 


0.7684 


4.5744 


100 


48.4 % 


0.2011 


1.8699 


1000 


56.5 % 


0.0557 


0.5327 


1 


% 


3.1088 


9.0965 


10 


12.3 % 


0.7775 


5.6462 


100 


54.4 % 


0.1620 


1.8172 


1000 


59.7 % 


0.0483 


0.5063 



TABLE I. Statistical properties of individual evolutions ac- 
cording to the Langevin scenario. The upper table belongs to 
the friction 77 = 2y p i , the middle to r\ — 7 P ; and the lower one to 
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